TTP10-12 

Narrow resonances studies with the radiative return methocQ 



O 

(N 

5h . 
Oh. 

< 

00 
(N 

S 1 

a; 

in: 

(N 
> 
OS 
t> 
(N 
O 

(N 
O 
O 



X 



Hcnryk Czyz, 1 Agnieszka Grzelihska, 2 and Johann H. Kiihn 3 

'institute of Physics, University of Silesia, PL-40007 Katowice, Poland. 
s Institute of Nuclear Physics Polish Academy of Sciences, PL-31342 Cracow, Poland 
s ' Institut fur Theoretische Teilchenphysik, Karlsruhe Institute of Technology, D-76128 Karlsruhe, Germany. 

(Dated: April 29, 2010) 

Using the radiative return method, experiments at high luminosity electron-positron colliders 
allow to explore the kaon and the pion form factors in the time-like region up to fairly high ener- 
gies. This opens the possibility to study kaon and pion pair production at and around the narrow 
resonances J/nj) and tp(2S) and explore the interference between electromagnetic and hadronic am- 
plitudes. Parameterizations of charged and neutral kaon as well as pion form factors are derived, 
which lead to an improved description of the data in the region of large invariant masses of the 
meson pair. These form factors are combined with the hadronic couplings of charged and neutral 
kaons to J/tp and ip(2S) and implemented into the Monte Carlo generator PHOKHARA, which is 
now, for the first time, able to simulate the production of narrow resonances and their decay into 
kaon, pion and muon pairs. 

PACS numbers: 13.66.Bc, 13.40.Gp, 13.20.Gd, 13.25. Gv 



I. INTRODUCTION 



New and precise measurements of the cross section 
for electron-positron annihilation into hadrons have been 
performed during the past years which were based on 
the method of " Radiative Return" [J, Q . Exclusive reac- 



tions, specifically two-body final states like tt+tt 



SI, 

pp [0|J or AA [6| and three- [7|, [8| and four-meson final 
states S [r| have been explored. An important ingredi- 
ent in these analyzes was and is the simulation of all 
these reactions through a Monte Carlo generator. In 
a first step, the generator EVA was developed 0, [Tl| . 
which is based on leading order matrix elements com- 
bined with structure function methods for an improved 
treatment of initial state radiation. Subsequently the 
complete next-to-leading order (NLO) QED corrections 
were evaluated [l2l . [l3l and implemented into the gen- 
erator PHOKHARA Tm4H, which is now available for 
a variety of exclusive final states. (For a recent review 
of theoretical and experimental results see e.g. [23|.) B- 
meson factories, operating at energies around 10 GeV 
and with high luminosity, allow to explore hadronic final 
states with relatively large invariant masses, up to 3 GeV 
and beyond. Therefore, the narrow resonances J/ip and 
^(2S) can be studied through the radiative return, in 
particular in decay channels of low multiplicity, leptonic 
ones like [241] . or two-body hadronic modes like 

tt+tt- , K+K~, K°K° or pp [|. The signal is identified 
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with the help of a very good mass resolution and particle 
identification in the resonance region. 



For an analysis exploiting the large statistics, the inclu- 
sion of radiative corrections from initial- and final-state 
radiation (ISR and FSR) is mandatory, since it affects 
the cross section and the line shape of the resonance. 
For the simulation of hadronic final states both the elec- 
tromagnetic contribution, i.e. a parameterization of the 
form factor, and the strength of the direct coupling of 
the resonance to the hadrons are required. The latter is 
absent for final states with positive G-parity (2tt, 4tt, ... 
) but non-vanishing e.g. for KK, 3tt or final states with 
baryons. On the other hand, a careful analysis of the res- 
onance line shape in the various channels would allow a 
model-independent determination of the direct coupling 
and of the form factors close to resonance I25rl3ll . 



With this motivation in mind we reanalyze the pion 
and kaon form factors with emphasis on the region above 
the p-resonance. The basic ingredients are very similar to 
those employed in an earlier study [32j ], However, addi- 
tional assumptions are required to properly describe the 
different resonance-like structures in the energy region 
between 1 GeV and 3 GeV. The details of this model 
and its parameters are described in sections II and III 
for pions and kaons, respectively. The new implementa- 
tion of these modes into PHOKHARA, which includes, 
as before, NLO ISR and FSR, is presented in Section IV. 
Section V is concerned with the implementation of the 
narrow resonances in the channels fi + \i~ , 7r + 7r~, K + K~ 
and K°K°. Hadronically and electromagnetically in- 
duced amplitudes are included, together with the radia- 
tive corrections from ISR and FSR. Section VI contains 
a brief summary and our conclusions. 
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II. 



THE PION FORM FACTOR 



For a realistic generation a model for the electromag- 
netic form factor is required. The ansatz presented in 

32} was published before the CLEO-c measurement of 
the form factor in the vicinity of the ip(2S) resonance 

33] and underestimates the experimental result signif- 
icantly. The same applies to the model predictions at 
JJjb as compared to the pion form factor calculated in 

31] from B{J/i)j — > tt + 7t~) and B(J/tp — > e + e~) decay 
rates. 

To accommodate the new data, the updated model 
ansatz for the pion form factor is taken similarly to [32j 



FAs) = 



N 



j fa 



n=(N+l) 



, (1) 



dQCD 



however, with different set of parameters. Those of the 
first N + 1 p radial excitations are fitted and the rest is 
taken from the "dual QCD model" [36j. It is necessary 
to take N = 5 to fit the data. For the precise treatment 
of P4 and p 5 see below. 

For the Breit-Wigner function we adopt the Gounaris- 
Sakurai [35j version with pion loop corrections included: 



BW Pn {s) 
where 
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s + H(s)-i^!T Pn {s) 



(2) 



H(s) = H(s) - H(m 2 J - (s - m 2 Pn )-H(m 2 pn ) , (3) 
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FIG. 1: (color online) The experimental data @, 0, OjJ [13- 
|42| | compared to the model fits results (see text for details). 
The form factor at J/ip comes from its theoretical extraction 
[3ll . 1 3 ll ] from the data. 



1 / 4m 2 

p(s) = -(s-4ml) 1 / 2 , v( s ) = ^/l-—?L. (5) 

Correspondingly we use the s- dependent widths 

which are taken from two-body P- wave final states and 
for simplicity (and lack of experimental information) also 



used for the rest of decay channels [32 ]. In Eqs. (jlj) and 
© we have used T Pn = T Pri (s = m 2 Pn ) , which is the total 
width of the p n meson. The constraint X)n=o° c p n = 1 
together with BW Pn (0) = 1 enforces the proper normal- 
ization of the form factor -^(0) = 1. 

For the ground state p(770) isospin violation from p—u 
mixing is taken into account by substituting 



3 



cl BW D0 (s) 

1 "T" C uj 

A Breit-Wigner function with constant width 

o 

= ^ — 



(7) 



(8) 



is used for description of the u resonance. 

As discussed in the Introduction, the couplings are 
based on the ansatz predicted in the dual-QC f Djv c =oo 
model 136] 



-i) n r(/3- i/2) 



V5fT(» + l)T(/3 - 1 - n) 



(9) 



where a' = 
tory a p (s) 



l/(2m%) 
-- 1 + a'(s 



is the slope of the Regge trajec- 
— mi ). The model postulates an 



equidistant mass spectrum mi n = m po (1 + 2n) and a lin- 
ear relation between mass and width of a given resonance 
r Pn = jm Pn , with 7 derived from the lowest resonance. 
The parameters (3 and m Po are to be taken from the fit. 

We fit the data in the time-like region which pro- 
vides detailed information about the structure of the res- 
onances and coincides with the region relevant for the 
PHOKHARA Monte Carlo generator. 

We have used new data |1, 0, |H, [H|-[42j , whenever 
possible. They are more accurate and the treatment of 
radiative corrections is well documented. Furthermore, 
we adopt the theoretical extraction of the pion form fac- 
tor at J/i/j using [3l[ 



AB(J/lp 7T+7T-) 




ftB(J/il> 



(10) 



and recent experimental data 



If one would assume independent point topoint statis- 
tical and systematic errors of the new data [3J, I38l442j and 
combine these in quadrature the results would be incon- 
sistent and no fit could be made. Summing linearly the 
statistical and systematic experimental errors for each 
experimental data point one finds very good agreement 
between the experimental data. This approach will be 
adopted below. The new BaBar data [J] become avail- 
able only after our analysis was finished and we include 
here only their part (above 1.2 GeV). The BaBar data 
below 1.2 GeV are in conflict with KLOE data and fur- 
ther investigations would be required how to merge these 
conflicting data samples. 

In [H, yH-liU the form factor including vacuum polar- 
ization was measured. We prefer to parameterize the 
'bare' form factor F v (see |25[ for definition), which is 
used throughout this paper and for example directly ob- 
tained in Eq. (|10p . The vacuum polarization corrections 



are taken from [20|, |43J . For the extraction of the form 
factor from the cross section, the CLEO-c collaboration 
'33} has corrected for the leptonic part of the vacuum po- 
larization effects. Hence their result has still to be cor- 
rected only for the hadronic part, which corresponds to 
a 1.5% shift of \F^\ 2 only and is irrelevant at the present 
experimental precision. 



Parameter 


model(fit) 


PDG value 
[34] 


model 


m P0 


773.37 ±0.19 


775.49 ±0.34 


input 


r 

1 PO 


147.1 ±1.0 


149.4 ±1.0 


input 




782.4 ±0.5 


782.41 ±0.12 


- 


r w 


8.33 ±0.27 


8.49 ±0.08 


- 


m pi 


1490 ±11 


1465 ± 25 


1340 


r« 


429 ± 27 


400 ± 60 


256 


m P2 


1870 ± 25 


1720 ±20 


1730 




357 ±46 


250 ± 100 


330 


m P3 


2120 [22] 


- 


2047 




300 [22] 


- 


391 


m Pi 


model 


- 


2321 


Tp 4 


model 


- 


444 




model 


- 


2567 


Tp 5 


model 


- 


491 


P 


2.148±0.003 




input 


K\ 


(18.7±0.5)-10 -4 






Arg{d) 


0.106 ±0.020 






|F 2 | 


0.59±0.10 






Arg(F 2 ) 


-2.20 ±0.16 








0.048 ±0.056 






Arg(F 3 ) 


-2. ±1.4 






\F 4 \ 


0.40 ±0.07 






Arg(F 4 ) 


-2.9 ±0.3 






\Fs\ 


0.43 ±0.05 






Arg(F 5 ) 


1.19±0.18 






X 2 /d-o.f. 


271/270 







TABLE I: Parameters of the pion form factor {Eg. ana 
Eq. Slip ) and results of the fit to the data. 



We have attempted to fit the experimental data keep- 
ing the coupling constants fixed to the model values 
(one fit parameter (3 for all of them) and fitting only 
the masses and the widths of the first few resonances 
(up to n = 5). This parameterization is satisfactory up 
to y/s ~ 1.3 — 1.4 GeV, where the details of the model 
for resonances, with n — 2 and higher, are not impor- 
tant. However, the model is definitively too simple for 
a description of the details of higher radial excitations 
including issues like coupled channels in decays of the 
higher radial p excitations. Hence we adopt a heuristic 
approach, where we allow for arbitrary complex couplings 
/„ of the p n (n = 1, 2, 3, 4, 5) 
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modulus square of the pion form factor and the form factor 
calculated with first six (N = 0, 5) resonances. 



(11) 



with F\ = 1 and four complex constants F n , n — 2, 3, 4, 5 
fitted to the experimental data. 



J 



The mass and the width of pz are fixed to their values 
obtained in the fit to the four pion production data [22j . 
For the masses and widths of the higher excitations (n > 
4) we use their model values. 

The results are shown in Table HI The fitted value 
of m po is smaller than its PDG2008 [3J| value,a conse- 
quence of using the dressed form factor in [3a - l42| . This 
phenomenon was also observed in pij . The parameters 
describing the radial p excitations obtained in the fit have 
to be taken with great care as they are strongly corre- 
lated, while in Table U we give only MINOS (MINUIT 
procedure from CERNLIB) parabolic errors. 

To ilustrate the numerical importance of the higher 
radial excitations within the "dual QCD model" in Fig. 
[5]we show the relative difference between the full modulus 
square of the pion form factor and the result calculated 
with the first six resonances which were used in the fit. 
It is evident that it is impossible to neglect the higher 
resonances and even in the po region they give small, but 
not negligible contribution to the form factor. 

III. THE KAON FORM FACTOR 



The kaon form factors were revisited for the same rea- 
sons as the pion form factor. Compared to the CLEO-c 
result [33| the model presented in [32[ underestimates 
the kaon form factor in the vicinity of the *p(2S) reso- 
nance. It is impossible to fit the existing data, inclu ding 
the CLEO-c result, with the functional form used in [32] 
or adding one or two more radial excitations, unless one 
would accept inclusion of a huge wide resonance in the 
region between J/ip and ifj(2S). To cure the situation, a 
model analogous to the one used for the pion form factor, 
assuming an infinite tower of resonances, was adopted. 
The ansatz reads 
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FIG. 3: (color online)The experimental data [33l. |45|. [4p-[5H compared to the model fits results (see text for details). The form 
factor at J/i/j comes from its theoretical extraction [2(| from the data. It was not used in the fit. 



The couplings in the part with subscript fit were fit- 
ted to the experimental data as well as the constants rj^ 
and c£. The values of N p ,N u and N$ are listed in Ta- 
ble [TTJ The entry PDG in Table [TT] implies, that masses 
and widths as given in PDG2008 [34[ were used. The 
masses and widths of the radial excitations, which were 
not measured, were calculated assuming an equidistant 
mass spectrum and a linear relation between the mass 
and the width of a given resonance 

m 2 ]n =rrij(l + 2n), T jn = 7 J -nij„ , j — p,uj,4> . (14) 

The value of 7 P was calculated from Eq.(fT4"f for n = 0, 
the other values were fitted to the data. 

Two versions of the model were investigated: the 'un- 
constrained' version were the couplings between kaons 
and p„, u> n and (j> n are not related and the 'constrained' 
version were = c¥ , n = 0, • • • oo. The 'constrained' 
model is not able to reproduce data as good as the 'un- 
constrained' model, however as evident from Tab. 1111 the 
corrections to the assumption c£ n = are small for the 
lowest two resonances. 

Despite this, the two models predict completely differ- 
ent asymptotic behaviour of the neutral kaon form factor 



in the region, where no data are available (Fig. @|. The 
'constrained' model, being closer to the SU(3) symmetric 
case where the neutral kaon form factor vanishes, arrives 
at significantly smaller predictions. The values of the 
couplings, which were not fitted, were calculated from 
the formula 

K _ {-impf - 1/2) 

° 0n ayim| n r(n + l)rOSf -1-n)' 

a' = l/(2m? ), j = p,u,(t> (15) 

with the exception of the couplings next to the last fitted, 
which were calculated from the normalization require- 
ments 

oo 

E c £ = 1 < 3 = P,u,4>- (16) 
n=0 

Breit-Wigner propagators 

9 

777 

BW c a = — ^— — , a = uj n (17) 

m£ — s — im a l a 
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with constant widths were used for all oj n , Breit-Wigner 
propagators with s-dependant widths 

m 2 , 

BWl 



ml - s - im<j, n T- 



3 ' 



r 



3 _ 



;i = K+,K° 



(18) 



were used for <fi n , the radial excitations of 0, and the 
GS Breit-Wigner functions (Eq.(j2|)) were used for p n . 
The parameters f3j , j = p,uj,(f>, were calculated from 
Eq.(fT5j) using fitted cj^ parameter. The results of the fits 
are summarized in Tab. HI] and in Figs. [3] and |U The 
high energy behaviour of both form factors is completely 
driven by the CLEO [33] measurement. 

Following |32j . i.e. assuming isospin symmetry, one 
arrives at the following predictions for the branching ratio 
of the r-lepton decay into K~K°i/ T 

Br(r~ — > K~ K a v T ) = (0.135/0.190) ■ (19) 



for 'unconstrained' and 'constrained' models, respec- 
tively. The model dependence is characterized by the 
spread between the two results and is far larger than the 
errors resulting from the fits of the parameters within one 
model. These results can be compared with the PDG 
value fH 



Br(r~ — > K K a v T ) = (0.158 ± 0.16) • 10" 3 (20) 
and are found to be reasonably consistent. 



Within the same assumptions one can predict the 
K~K° invariant mass distribution and compare it (see 
Fig. [5]) with existing CLEO data [Hj]. As evident from 
Fig. [5]both models give very similar predictions and both 
agree with the data. Thus we conclude that within the 
current experimental accuracy, which is, however, very 
poor, isospin symmetry works well and the details of the 
models do not play any role in its tests. 



IV. MONTE CARLO IMPLEMENTATION OF 

K+K~ AND K°K° 

The event generator PHOKHARA has been extended 
to generate K + K~ and K°K° final states. In this section 



we present the implementation and results for the region 
below the narrow resonances J/tjj and ip(2S). Charged 
kaons have been implemented in the same way as the 
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Parameter 


Input 


Fit(l) 


Fit(2) 


PDG value 


model(l) 


modcl(2) 


m <Po 
r 0o 


: 


1019.415 ±0.004 
4.34 ±0.01 


1019.415 ±0.003 
4.22 ±0.04 


1019.455 ±0.020 
4.26 ±0.05 


input 
input 


input 
input 


m^j 
r 0i 


1680 
150 






1680 ± 20 
150 ±50 


1766 
353 


1766 
353 


m P0 


775.49 
149.4 


; 
- 


; 
- 


775.49 ±0.34 
149.4 ±1.0 


input 
input 


input 
input 


m pl 

r„ 


1465 
400 






1465 ±25 
400 ± 60 


1345 
259 


1345 
259 


r 

1 P2 




1680 ±4 
365 ± 59 


PDG 
PDG 


1720 ± 20 
250 ± 100 


1734 
334 


1734 
334 




782.65 
8.49 


- 


- 


782.65 ±0.12 
8.49 ±0.08 


input 
input 


input 
input 


r 


1425 


145 ±9 


PDG 


1400-1450 
180-250 


1356 
678 


1356 
678 


771^2 




1729 ± 76 

245 ±9 


PDG 
PDG 


1670 ± 30 
315 ±35 


1750 
875 


1750 
875 






1.040 ±0.007 


1.055 ±0.010 








70 

r K 

00 
„K 
C 01 

„K 

C 02 

„K 

C 03 


- 


1.97 ±0.02 
0.2 

0.985 ±0.006 
0.0042 ±0.0015 
0.0039 ± 0.0026 
(JTU)0.0033 ± 0.0067 


1.91 ±0.02 
0.2 

0.947 ± 0.009 
0.0136 ±0.0024 
(00.0214 ± 0093 


- 
- 


input 
input 
0.0084 
0.0026 
0.0012 


input 
input 
0.0271 
0.0088 


OO 

E 4 n 


model 


0.0036 


0.0180 




0.0036 


0.0180 


Hp 

7p 
U P1 

P2 

c K 

P3 
K 


en 

- 


2.23 ±0.06 
0.193 dUll (= T p /m p ) 
1.138±0.011 
-0.043 ±0.014 
-0.144 ±0.015 
-0.004 ±0.007 
(JTU)0.0662 ± 0.0243 


2.21 ±0.05 
0.193 JUJ) (= r p /m p ) 
1.120 ±0.007 
-0.107±0.010 
-0.028 ±0.012 
(|j"6l)0.032 ± 0.017 


- 
- 


input 
input 
-0.087 
-0.020 
-0.0084 
-0.0045 


input 
input 
-0.078 
-0.019 
-0.0079 


OO 

E <& 

n = JV p + l 


model 


-0.0132 


-0.0170 




-0.0132 


-0.0170 


7" 

c5 
t& 

OO 

E <4 

n=JV„ + l 


c£ mi 
model 


Pp 
0.5 

c K 

PO 

c K 

P2 

c K 

PS. 

L P4 

oc 

E <ft 

n = N p + l 


2.75 ±0.06 
0.5 

1.37 ±0.03 
-0.173 ±0.003 
-0.621 ±0.020 
([T6ll0.43zfc0.04 

-0.0096 




input 
input 
-0.087 
-0.020 
-0.0084 
-0.0045 
-0.0132 


input 
input 
-0.345 
-0.026 
-0.0079 

-0.0096 


X 2 /d-o,f. 




277/256 


221/260 









TABLE II: Parameters of the kaon form factors and results of the fit to the data. Masses and widths are given in MeV. The 
column 'Fit(l)' (Fit(2)) contains the values of the constrained (unconstrained) fits. 
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FIG. 5: (color online) 

dT (t -> K ~ it i/ T ) / dJ~Q 



Normalized distributions 
of the kaon pair invariant mass pre- 



r(r->K- K°v T ) 

dieted within two models, described in the text. CLEO data 
55] normalized to the total number of events are also shown. 



7r + 7r~ channel [HI, with the kaon form factor de- 
scribed in Section Mil The NLO FSR corrections have 
been implemented as well. For the neutral kaons the 
corrections are limited to ISR. With the enormous lu- 
minosity of B factories one expects hundreds of events 
even for Q 2 between 3 and 4 GeV 2 and large statistic 
around the resonance (Fig. [SJ. The next-to-leading 
FSR corrections are relevant for a measurement in the 
neighbourhood of the <f> resonance if an accuracy better 
then 10% is aimed (Fig. [7}. 
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FIG. 6: (color online) Differential cross section for y/s = 
10.52GeV of the processes e + e~ — > K + K~ 7(7) and e + e~ — > 
K° K°j( / y) with angular cuts. 



V. NARROW RESONANCES AND THE 
RADIATIVE RETURN 

The radiative return, as implemented now in 
PHOKHARA, receives contributions from a multitude 
of amplitudes shown in Fig. [5J The notation introduced 
in this figure also applies to the narrow resonance am- 
plitudes. Below we list the formulae for 7r + 7r~, 
K+K- and K°K°. For the K°K° amplitude FSR emis- 
sion is not present. We explicitly indicate the vacuum 
polarization contributions as well as the <f> contribution. 
We assume that in the vacuum polarization J/ip, ip(2S) 
and <f> were not accounted for. The differential cross sec- 
tion given by PHOKHARA reads 



da = 

\M 1uLO isr ■ Cl P P (Q 2 ) + M- (1 ,lofsr ■ C^ p (s)\ 2 d^ 
+ \M 2lJ sR-Cl P P {Q 2 )\ 2 d<S> 2 

+2 MM 1u nloisr x Ml i LOISR ) ■ \C^ p P (Q 2 )\ 2 d^ 

+\M luISR , ^ FSR ■ Cl P P {{Q + fc 7 ) 2 )| 2 d$ 2 

+2 Re(M^ L %r R x M\ uLOISR ) ■ \C v R p P {Q 2 )\ 2 d^ 

+ \M lJS R , lu FSR ■ Cl P p ((Q + fc 7l ) 2 )| 2 d$ 2 

+2 Re(M»%™X R x M\ uLOFSR ) ■ \C^ p P (s)\ 2 d^ , 

(21) 

where 
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FIG. 7: (color online) Comparison of IFSNLO and ISRNLO 
distributions with angular cuts. 
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Ml 



(22) 



(23) 



and (d$2) denote the phase space with one (two) 
photon(s) in the final state with all statistical factors 
included. 

For P = fj, and P = ir, Cp = (no direct decay of the 
narrow resonances into n + fJ-~ and tt + tt~), while 5p = 
for P — K and 5p = 1 for P = /j, and P = ir. The 4> con- 
tributions to the kaon pair production are included in the 
kaon form factor, hence 8k = 0. The notation and the 
detailed description of the narrow resonance contribution 
to the amplitude can be found in (2{| (see [13] for simi- 
lar studies). From [|| we also take \c J Jf\ = 1.27 ± 0.32 

and \c^l S) \ = 2.94 ± 0.99. The information on the neu- 
tral kaon couplings to the narrow resonances is almost 
nonexisting and we use the lower limits of |c^'| = 2.81, 
| c ^(2S)| _ g which correspond to the upper limit on 
the neutral kaon form factor (see [25[ for details). The 
phases are essentially not known and we use 100° to ob- 
tain the numerical values in the next section. 



VI. THE IMPLEMENTATION OF NARROW 
RESONANCES INTO THE MONTE CARLO 
EVENT GENERATOR PHOKHARA 

The tests of the ISR part of the implementation of the 
narrow resonances where straightforward and followed 





\ iV1 jiJSR ; -y,FSR\ 



.WISH' 





\MjjSR ; ji,F3r\ 



FIG. 8: (color online). The contributions to the radiative 
return cross section included in PHOKHARA. Label '1' at a 
photon line means that the photon is hard. 



the standard tests we perform for each new channel 0- 
|20| . HH . The comparisons were made with the analytic 
formulae of [53j separately for one and two photon emis- 
sion. The precision of the comparisons was at the level of 
a small fraction of a per mill, proving the technical pre- 
cision of the program at that level. The independence of 
the results on the separation parameter between soft pho- 
ton, calculated analytically and hard photon, generated 
by means of the Monte Carlo method, was also tested 
with that precision. 

The implementation of the NLO FS R p art is more 
tricky. The analytic formulae used in [l6|, Il8| for soft pho- 
ton contributions are still valid, which we have checked 
numerically with the precision of 0.02%. However if one 
chooses the separation parameter between soft and hard 
part at the usual value w = 10 ~ 4 , which corresponds 
to the photon energy E~ t = 1 MeV for y/s = 10 GeV, 
the 'soft' integral receives contributions from the whole 
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FIG. 9: (color online) Comparison between invariant mass 
distributions obtained with w = 10~ 4 and w = 10~ 7 . 
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FIG. 10: (color online) Comparison between invariant mass 
distributions obtained with w = 10~ 4 and w = 1CF 7 smeared 
with Gaussian with the standard deviation of 14.5 MeV. 



resonance region, as a consequence of the small width ( 
r = 93.4 keV). For a cutoff of 10~ 4 the part of the ma- 
trix element, which multiplies the soft emission factor is 
rapidly varying and the basic assumption underlying the 
whole approach, that the soft emission can be integrated 
analytically with the multiplicative remainder being con- 
stant, is not longer valid. Pushing the value of the cutoff 
to an extremely small value, say 10 -7 , solves this prob- 
lem. However, single-photon emission is not an adequate 
description for such soft photons and in principle one 
should use exponentiation. From the technical side this 
is reflected in the appearance of negative weights. In- 
clusion of YFS-like multi-photon production would allow 
to cure this problem. However since this would amount 
to completely restructuring our Monte Carlo generator, 



FIG. 11: (color online) Comparisons of the invariant mass dis- 
tributions obtained with w — 10 -7 taking into account only 
ISRNLO contributions and the complete (ISR+FSR)NLO re- 
sult. 
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FIG. 12: (color online) Comparisons of the invariant mass dis- 
tributions obtained with w = 10~ 7 taking into account only 
ISRNLO contributions and the complete (ISR+FSR)NLO re- 
sult. Detector smearing effect are taken into account. 



we have adopted a simpler approach, which gives correct 
distributions, when convoluted with an energy resolution 
typical for a detector at a <fi- or B-meson factory. 

Due to the finite detector resolution one never observes 
the true distribution of the events, but the one convoluted 
with the detector resolution function. This increase of 
the effective width by about a factor hundred is sufficient 
to cure the problem. For a cutoff of 10~ 4 the distribu- 
tion remains smooth and we can produce the unweighted 
events sample. The result will, as expected, depend on 
the resolution of the detector. To check if this is true for 
the realistic energy resolution of the BaBar detector [24| 
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FIG. 13: (color online) Relative ratio of the invariant mass 
distributions taking into account only ISRNLO contributions 
and the complete (ISR+FSR)NLO result. Detector smearing 
effect are taken into account. 



of 14.5 MeV we have compared the muon invariant mass 
distributions obtained with w = 10~ 4 and w — 10~ 7 , 
smeared with a Gaussian distribution with a standard 
deviation of 14.5 MeV. Even if the non-smeared distri- 
butions are completely different, as shown in Fig. [5] the 
smeared distributions agree within 2 per mill as shown in 
Fig. 1101 This 2 per mill is the intrinsic error coming from 
the method we use, but the generator should be accurate 
enough for any practical purposes. 

It is interesting to observe (Fig. QT]) that the FSRNLO 
contributions fill completely the interference dip, still vis- 
ible if only ISR corrections are taken into account. Thus 
the absence of the dip in the observed invariant mass dis- 



tribution is not only the effect of the detector smearing. 

The huge FSRNLO corrections seen in Fig. [TT] arc 
washed out if one looks at the detector smeared distri- 
butions shown in Fig. 1121 The corrections are seen more 
accurately in Fig. 1131 where the relative difference is 
shown. The FSRNLO corrections cannot be neglected if 
one aims at a precision better then 10%, unless one con- 
siders only the integral over the whole resonance region 
(together with the side bands as shown in Fig. [T"2"|) . In 
the integrated cross section a large part of the corrections 
cancel {visrnlo = 6.901 pb, uifsrnlo = 6.954 pb). 

Identical tests were performed for the pion pair pro- 
duction with identical conclusions, so we do not present 
them here. 



VII. 



SUMMARY 



New parametrizations of the pion and kaon form fac- 
tors, based on the "dual QCD model", are presented 
which are derived from a fit to a combination of old 
measurments and more recent experimental results in the 
energy region above the p-resonance. These form factors 
and the results of a recent analysis of the direct hadronic 
coupling of KK to J/tp and i){2S) are incorporated in 
the Monte Carlo generator PHOKHARA which is now 
also adapted to the simulation of narrow resonances, in- 
cluding the effects of ISR and FSR in NLO. 
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